UWThPh-2006-1 



An open-quantum-system formulation 
of particle decay 

Reinhold A. Bertlmann,* Walter Grimus,^ and Beatrix C. Hiesmayr 

Institut fur Theoretische Physik, Universitat Wien 
Boltzmanngasse 5, A-1090 Wien, Austria 

February 14, 2006 



Abstract 

We consider an open quantum system which contains unstable states. The time 
evolution of the system can be described by an effective non-hermitian Hamiltonian 
H e g, in accord with the Wigner-Weisskopf approximation, and an additional term 
of the Lindblad form, the socalled dissipator. We show that, after enlarging the 
original Hilbert space by states which represent the decay products of the unstable 
states, the non-hermitian part of H e Q — the "particle decay" — can be incorporated 
into the dissipator of the enlarged space via a specific Lindblad operator. Thus 
the new formulation of the time evolution on the enlarged space has a hermitian 
Hamiltonian and is probability conserving. The equivalence of the new formulation 
with the original one demonstrates that the time evolution which is governed by a 
non-hermitian Hamiltonian and a dissipator of the Lindblad form is nevertheless 
completely positive, just as systems with hermitian Hamiltonians. 
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Introduction: In reality a quantum system is not isolated but always interacting with 
its environment and has to be considered as an open quantum system [H 12] • It leads to 
a mixing of the states in the system — decoherence — and to an energy exchange between 
system and environment — dissipation jSl IH 13 Ej • The decoherence/dissipation weakens 
or destroys the typical quantum phenomena, the interferences. 

There is great interest to study such decoherence/dissipation processes for elementary 
particles like the "strange" K-mesons (neutral kaons) and the "beauty" B- mesons 0011]], 
in particular, the decoherence of entangled meson pairs [TTH ITT) IT2*] (for an overview 
see Ref. ^S])- However, these particles are decaying, which leads in the framework 
of quantum mechanics to some difficulties. In the usual framework the decaying states 
evolve according to the Wigner-Weisskopf approximation (WWA) where the probability 
of detecting the particle is not conserved and the corresponding Hamiltonian is, for this 
reason, non-hermitian. 

The master equation for such an open quantum system described by the density matrix 
p{t) is given by 

^ = -iH cS p + ipHl s - D[p], (1) 
where the dissipator has the following general structure of operators [Tl) IT5] : 

DIP] = \i: (A)A 3 p + pA)A 3 - 2A jP A\) . (2) 

j 

The operators Aj are usually called Lindblad operators. The effective Hamiltonian of the 
system is given, according to the WWA, by 

H cS = H- l -r, (3) 

where H and T are both hermitian and, in addition, T > 0. 

In terms of open quantum systems, the time evolution of the density matrix represents 
a dynamical map which transform initial density matrices p(0) on the Hilbert space of 
states, H s , to final density matrices pit) = V t [p{0)], while the system is interacting with 
an environment. Such dynamical maps are i) trace conserving, ii) convex linear, iii) 
completely positive. Complete positivity is a rather strong and important property. It 
is defined by demanding that all extensions V t ® t n on the Hilbert space H s ® <C™ are 
positive, i.e. 

(T4®l„)[p]>0 Vt>0, Vn = 0,l,2,..., and Vp>0, (4) 

where p is a density matrix on H s <g) <C n . Physically, it is a reasonable condition since 
the extension Vt ® l n can be considered in quantum information as an operator on the 
composite quantum system Alice and Bob, which acts locally on Alice's system without 
influencing Bob. Furthermore, complete positivity essentially ensures that tensor products 
of maps Vt remain positive, an important property especially when considering entangled 
states. 

However, extracting from a dynamical map the Hamiltonian term in the usual way 
leads to a hermitian Hamiltonian in a master equation like Eq. (Q). In this sense the 
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authors of Ref. ^7j investigated recently decoherence and decay of the neutral kaons by 
using a dynamical map corresponding to a master equation with hermitian Hamiltonian. 
Nevertheless, many authors work with a non-hermitian Hamiltonian — which is a standard 
procedure in particle physics (see, e.g., Ref. (THj) — in order to include the decay property 
of the particle. 

The formalism: We are going to show in this letter that both approaches — each one 
having its own appeal — are indeed equivalent to each other for quite general quantum 
systems. To be more precise, we ask the following question: 

Can we work with a hermitian Hamiltonian and incorporate the decay as Lindblad operator 
such that the time evolution of the system represents a completely positive map? Moreover, 
does it describe the WW A properly without effecting decoherence and/ or dissipation? 
The answer will be yes, it is possible! But we have to enlarge the Hilbert space and 
include formally the decay products of the unstable states. 

Thus we assume that the total Hilbert space H to t is the direct sum 

H tot = H s ©H / , (5) 

where H s contains the states of the system we are interested in whereas Hf is the space 
of the "decay states," defined in the following way. We assume that d s = dimH s < oo. 
Then the non-hermitian part T of the effective Hamiltonian in Eq. (jH)) can be decomposed 
as 

r 

r = I^Tj l^'X^'l w^h 7 i >0Vj, (6) 

3=1 

with an orthonormal system {fj} in H s and r = dimH s — n , where no denotes the 
degeneracy of the eigenvalue zero of T. We demand that dimHj = dj > r, otherwise Hy 
is arbitrary. 

In the following, p denotes a density matrix on H tot , where it has the following de- 
composition: 

9 = ( P'fs p'ff ) P " = PSS ' p]ff = P// ' = PSI ' 

Now we have to define the time evolution p on H tot . The Hamiltonian H and the Lindblad 
operators Aj are easily extended to the total Hilbert space by 

»-(?!!)■ (8) 

Note that for the definition of 7i we have used the hermitian part H of H c g. 

Now we turn to the decay. We need a Lindblad operator B on the full space H tot 
which describes the decay in the subspace H s . As we will see in the next paragraph, the 
decay is described by 

8 = { r n ) with B ■ H ^ H / ( 9 ) 
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and 

Y = B ] B. (10) 
Let {fk} be an orthonormal basis of H^. Then we can decompose B as 



a = EE 6 wl/*>tel- ( n ) 
fc=i i=i 

In order to fulfill Eq. we require 

«*/ 

J2 b *kihj = Sijjj- (12) 
fe=i 

For df > r, such a c?/ x r matrix (6^-) always exists. The simplest case is df = r, where 
we can choose 

df 

£ = EVW;X^I- (13) 

In that case, each unstable decaying state \(pj) would decay into just one specific decay 
state \fj). 

Then we can write the following master equation for the density matrix on H tot : 

^ = -i[n,p] -v[ P ], (i4) 

with the dissipator 

V[p\ = ^ E (M'P + P A \ A i ~ 2 AM]) + \ { B]B P + P& B ~ 2BpB f ) • (15) 
2 j 2 

To prove that this master equation contains Eq. (JTJ), we decompose it with respect to the 
components of p: 



p s/ = -iH Ptf - -B^Bp a f - -^A^AjPsf, (17) 
p,/ = Sp ss 5t. (18) 



Indeed, with Eq. (110)1 . we immediately see that Eq. ()16)) for p ss reproduces the original 
equation (0). Furthermore, it is obvious from Eqs. (fT4*j) and (|T5^) that Tr p(t) = 1 Vt > . 

At this point, some general remarks are at order. By construction, the time evolution 
of p ss is independent of p s f, Pf s , and pff. Actually, the time evolution of p s f or pf s 
completely decouples from that of p ss and pff. The time evolution of p//, the density 
matrix of the decay states, is determined solely by p ss {t) and is the characteristic for a 
decay process! 
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With the initial condition p//(0) = 0, it is simply given by 

p ff (t) = B f&t' Pss {t')B\ (19) 

J 

If we choose the initial condition p s /(0) = 0, then p s f remains zero for all times and the 
same applies to pf s . Anyway, these parts of p are totally irrelevant for our discussion. 

Generally, for decaying systems the following properties of the time evolution (fPfj) are 
most important: 

• p ss (t) and Pff(t) are positive Wt > 0; 

• The time evolution of p ss is completely positive. 

Proof: The density matrix p(t) obeys the time evolution given by the master equa- 
tion (JHJ), therefore, p{t) > Vt > 0. Thus, confining ourselves to vectors v G H s , we 
evidently have < (v\p(t)v) = {v\p ss (t)v) and, therefore, p ss (t) > 0. The same reasoning 
holds for Pff(t). As for complete positivity, we note that on H tot ®C a general operator 
has the structure R = J2k ^fc^-^fc; with operating on H tot and L k on C". We know that 
the dynamical map V t induced by Eq. (|14j) is completely positive. Therefore, if R > 0, 
then we have (V t ® t n )[R] > 0. We also know from Eqs. (fTKll. (fTTfl . and (JHJ), that V t has 
a well defined restriction to operators acting purely on H s . Thus, R > on H s ® C n is 
mapped into a positive operator on the same space by V t <g> l n . This concludes the proof. 

The case of a non-singular T: Finally, it is interesting to consider the special case of 
a non-singular T, where all states of H s decay. Then, dimH^ > dimH s and, in addition 
to what we discussed before, the following properties of the time evolution (fTlj) hold: 

i) lim t ^ 00 p s/ (t) = 0, 

ii) lim^^ Pss (t) = 0, 

iii) lim^ooTrp//^) = 1. 
Proof: 

i) We consider the equation 

— uj = -iHu - -V'u, with Y' = Y + YA\A j (20) 
dt 2 j 

and cu(t) G H s . From Eq. (J2"U|) we derive 

A|^|2 = _ (w |ra;)<- 7 ;|a;| 2 , (21) 

where 7q > is the smallest eigenvalue of T'. Since \u\ 2 > 0, it follows from this 
equation that |w(t)| 2 < |c<j(0)| 2 exp(— 7g£) and, consequently, lim^ 00 a;(t) = 0. This 
leads to lim^ 00 p s /(t) = 0, because p s f is a linear combination of operators 
with uj G H s , / G H/. 
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ii) Taking the trace of Eq. we obtain 

^Trp ss (t) = -Tr (Tp ss (t)) < - 7o Tr p ss (t) . (22) 

Here, 70 is the smallest eigenvalue of T. From Eq. (}2"2"|) we derive the inequality 
Tr/9 SS (£) < Tr p ss (0) exp(— 7 t), whence we find lim^oo Tr p ss (t) = 0. Since p ss (t) > 
0, the whole matrix must vanish in the limit t — > 00. 

iii) Point ii) and Trp(t) = Tr p ss (t) + Tr pff(t) = 1 imply lim^oo Tr Pff(t) = 1. 



The simplest possible example: We choose the minimal dimensions d s — df — 1, 
the Hamiltonian 

«=(?° ). <») 

and set Aj = 0, i.e. we neglect any decoherence part. Then, Eqs. (fTBj). (|TH|) and (fT^j) 
simplify considerably, providing the result 




where we have chosen the initial conditions p ss (0) = 1, p//(0) = p s /(0) = 0. This example 
illustrates the general fact that, for non-singluar T, the probability loss of the system is 
balanced by the probability increase of the decay states. 

Considering the mixedness of the quantum states, we can use as measure 

S(t) = Trp 2 (t) = l-2e- Tt + 2e- 2Tt . (25) 

We have 6(0) = 1, lim^oo 5(t) = 1 and 5(t) < 1 for < t < 00, i.e. in the beginning, p 
represents a pure state, then the state is mixed, whereas for large times it approaches a 
pure state again, when the system definitely has changed into the decay state. 

The case of neutral kaons, where the Hamitonian H and the operator B are 2x2 
matrices, has been investigated in detail in Ref. [T7| . 



Kraus operators: Finally, in the case of our example, we would like to consider the 
time evolution of the density matrix as a dynamical map V t . In our case, V t can be 
represented by a sum of operators [HI2] : 

p(0) — V t \p(0)] = £ Mj(t) p(0)Mj(t) = p(t), (26) 

3 

with the normalization J2j Mj(t)Mj(t) = 1 for the Kraus operators Mj(t). For the corre- 
spondence between the Lindblad and Kraus operators see Ref. [TH] . 

The operator sum representation (j2T?j) is a very useful approach in quantum information 
to describe the quantum operations or the specific quantum channels. In our example, 
the decay of a particle corresponds to the amplitude damping channel of a quantum 
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operation [THj, e.g., the spontaneous emission of a photon, and we can calculate the 
Kraus operators needed: 



with the probability p(t) = 1 — e~ rt . 

Conclusions: In this letter we have considered a time evolution given by Eq. (0), with a 
non-hermitian Hamiltonian H e R and a dissipator of the Lindblad form (J2J). Assuming that 
the non-hermitian part of H e & describes "particle decay," we have shown that such a time 
evolution is completely positive. Our strategy was to add the space of "decay states" H/ 
to the space of states H s and to extend the time evolution in a straightforward way to the 
full space H s © H/, such that this full time evolution is probability-conserving. With the 
initial condition that the density matrix is only non-zero on H s , the Lindblad operator B 
of Eq. ©, which is responsible for the decay, shifts states from H s to Hj, whereas at the 
same time the Aj terms cause decoherence and dissipation on H s . Thus, particle decay 
and decoherence/dissipation are related phenomena, which can be described by the same 
formalism of a completely positive time evolution. 
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